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It has been known for a long time that perturbative calculations in hot SU(N) gauge 
theories(QCD) are plagued with an infrared problem. The scale of momentum where the 
problem emerges is g 2 T where g is the gauge coupling and T is the temperature of the 
plasma. No reliable way has been suggested for treating these long range modes. The 
problem can be seen as originating from the fact that thermal fluctuations of the magnetic 
modes at the scale g 2 T are large enough for the theory to become essentially non-linear. 
The g 2 T scale is important for understanding a large number of physical problems such 
as magnetic screening!] baryon number transitions in the electroweak theory at high 
temperature and the determination of the parameters of the electroweak phase transition. 
In particular, all of these need to be understood for a complete calculation of electroweak 
baryogenesis. The real-time treatment is appropriate for some problems, especially for 
computing topological transition rates in the Standard Model. 

In this paper we derive a set of equations which describe the real-time dynamics of 
non-perturbative soft modes in hot gauge theories and we discuss two of their most im- 
portant applications. The first one is baryon number violation at high T. The physics 
of the soft modes is stochastic in nature and has the natural time scale ~ (g^T)^ 1 . In 
particular, our equations suggest that baryon number violation at high temperature in 
the unbroken phase scales as ajj^T 4 as the result of a diffusive motion across a topolog- 
ical barrier, rather than a^T 4 as classical arguments and recent numerical simulations 
suggest 0, [|. Last but not least, our equations constitute a solid computational tool 
for the determination of the numerical coefficient of a^T 4 . This would have important 
consequences for electroweak baryogenesis, and the constraints imposed on it from the 
structure of the symmetry breaking sector and from the CP violating sector. As a second 
application, we introduce a possible path for the real-time resolution of infrared diver- 
gences which plague the perturbative computation of the parameters of the electroweak 
phase transition. 

Stochastic dynamics of soft modes in hot gauge theory 

1. Let us begin with a simple argument showing that the time scale responsible for 
processes occurring at the space extent of {g 2 T)~ 1 has the natural time scale of (g 4 ! 1 ) -1 .^] 
To this end we consider thermal fluctuations of the gauge field A^. The power spectrum 
of these fluctuations can be related to the imaginary part of the gluon propagator by the 
fluctuation-dissipation theorem. At low frequencies, w CT, one has, 

2T 

(A;(u, q) q')> = hnD*{u>, q) • (2n) 4 S(co - <S)6{q - q') (1) 

to p 

1 For recent attempts to compute the magnetic mass, see fl and references therein. 

2 For parallel considerations with an emphasis on cross barrier transitions see . 
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II(u;,q) = -icg 2 T 2 



The retarded propagator D R at small u and q acquires important contributions from the 
hard thermal loops, 

D R (u, q) — r . 

uj l — q^ — ll(u;, q) 

For non-perturbative physics, we are interested only in transverse (or magnetic) modes, 
for which the gluon polarization II(a>, q) has the following behavior 

q 

in the regime u <C q = |q|, where c is some real constant. Eq. ([I]) implies, in this regime, 

<A>, q)A{u/, q')> ~ g6 /^y^ • ^ " - q') • (2) 

To compute the one-time correlation of the field, we integrate Eq. (g) with respect to u>; 
this gives 

(A*(t, q)A(t, q')> ~ 5(q - q') / du (-l)-Im- i , ~ l*(q - q') (3) 

where the integral is saturated for uj ~ g~ 2 T~ 2 q 3 . Substituting q ~ g 2 T in this expression 
implies uo ~ g 4 T, which is what we wanted to show. To see that the non-perturbative 
scale of spatial momenta is g 2 T, we can use Eq. (^) to find the contribution of the modes 
with momenta of order q to the correlator of A at coinciding points; we find 

(A 2 (t,x))~Tq. 

Hence the typical fluctuation size of A is (Tg) 1 / 2 if restricted to modes with spatial momen- 
tum of order q. The nonlinear contribution to the field tensor gA 2 becomes comparable 
with the linear part dA when q ~ g 2 T. So, we have shown that the non-perturbative 
physics is associated with modes with spatial momenta of order g 2 T and frequencies of 
order g 4 T. 

2. To find the effective theory describing these non-perturbative modes, we resort to 
the recent reformulation of the hard thermal loop Lagrangian in terms of a kinetic Vlasov 
equation [|5], ^J. In QCD, the soft dynamics is described by the soft field -4^(2, x) and 
8n a (t, x, p), the deviation from equilibrium of the hard gluon distribution functions. The 
scale of space and time variation of both A^t, x) and Sn a (t, x, p) is much larger than T -1 
while p ~ T. The time evolution of these two effective variables of the theory is described 
by a Langevin-type Vlasov equation having the form 

v»D*5n% x, p) + gv • E a ^|^ = (4) 
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with 

d»F; u = 2gN J vjn a (t, x, p) + Q(t, x) (5) 

where w M = (l,p/|p|), and D a ^ = 8 + gf abc A c is the covariant derivative. In these 
equations, n refers to the equilibrium Bose-Einstein distribution. The difference between 
Eq. (|5|) and that discussed in Ref. || is the presence of the stochastic noise £"(x) in 
the field equation. Eqs. could be derived using the formalism of Ref. |7j which 

has recently been implemented in the scalar theory ||. However, in this letter we will 
rather use physical arguments to justify the need to include the stochastic noise into 
the kinetic equation and to find the correlator (Cu( x )^v(v)) ■ m f ac ^ while the noiseless 
Vlasov equation describes the time evolution of the main values of the mean field A^(x) 
and the distribution function <5n a (t,x, p), there are many instances where spontaneous 
thermal fluctuations need to be considered (for example, in the problem of baryon number 
violation in a hot plasma). Because the term in Eq. (]5|) reflects fluctuations in the 

current density, one of our main tasks is to find the form of the noise correlator and to 
simplify Eqs. (flJ^) so that they include only the modes we are interested in, i.e., modes 
with u} ~ g 4 T and p ~ g 2 T. 

Before proceeding further, let us make one important remark. Eqs. (f|,|5|) have been 
shown in Ref. || to reproduce the HTL (hard thermal loop) effective Lagrangian describ- 
ing the physics of modes with momentum q ~ gT. One could ask whether the HTL 
Lagrangian correctly describes the physics at much smaller spatial momentum scale of 
g 2 T. A diagrammatic analysis shows that beside the HTL diagrams, diagrams with ladder 



insertions are also important (see and a recent discussion in |10|). However, we argue 
here that these diagrams do not modify the dynamics at the momentum scale g 2 T. Our 
first argument is based on the similarity of the ladder diagrams with those considered in 



TT| in the framework of the scalar theory. In the latter context, it has been shown that 



the effect of this infinite set of diagrams is the generation of a collision term in the trans- 
port equations. One can argue that this statement can be extended also to gauge theories 
[ I2| . Therefore, we can take into account all diagrams relevant to the g 2 T dynamics by 
simply adding a collision term to the RHS of Eq. (§). However, as will be demonstrated 
later, the effect of collisions is suppressed at length scales (g 2 T)~ x and starts to becoming 
important only at length scales of order (g 4 T) _1 . As a bonus, we can safely ignore all 
non-HTL diagrams. A second argument is based on explicit calculations in QED show- 
ing that the resumed diagrams do coincide with the naive HTL answer for U(u>,q) when 

3. Let us now derive the correlator of We begin by considering in detail a 
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simpler case, namely, hot QED, for which the Vlasov equation has the form, 



g + v |) Mi , x , p) + eEv MM.o (a) 

PF^ = 4e J ^-v„5n(t, x, p) + £,(*, x) . (7) 

Here 8n(t, x, p) is the deviation from equilibrium of the distribution function of fermions 
with a given spin (we assume that the distribution functions are diagonal in the spin 
space), and n is the equilibrium Fermi distribution. The factor 4 is due to the fact that 
the plasma consists of particles and anti-particles, each having 2 spin degrees of freedom.^] 
The dynamics in QED is simplified, since the Vlasov equation is linear, but still, it is 
useful to consider it before turning to the more complex case of QCD. 

To see the physical origin of the stochastic noise £„ in Eq. (|7j), let us imagine the plasma 
as a collection of fermions and anti- fermions, with, for distribution functions, n + (t,x, p) 
and n_(i, x, p), respectively. In the absence of an external field, n±(t, x, p) are equal to 
the thermal distribution n(p) and the resulting current density in the plasma vanishes 

j,(t, x) = e J-^- 3 v, J>+(*. x, p) - n1(t, x, p)) = , (8) 

(the index a labels the spin degrees of freedom). However, Eq. (H) represents only the 
average of the current density: there are always thermal fluctuations of j^(t,x.) due to 
fluctuations of the distribution functions n±(t, x, p) about the Bose-Einstein main value. 
These fluctuations can be divided into a long-distance part and a short- distance part. The 
former is 5n(t, x, p) which enters the Vlasov equation (Eq. (^)). These fluctuations have 
for typical length scale ~ {g 2 T)~ 1 . The short-distance part fluctuates over much shorter 
distance scales (up to T _1 ). These short-distance fluctuations are essential for keeping 
the soft modes in thermal equilibrium. A Vlasov equation without the stochastic noise 
predicts that modes with u < q go away after a sufficiently large amount of time due to 
Landau damping: this is not the case of fluctuations in a realistic plasma. Let us denote 
the short-distance fluctuations of the distribution function by An a (t, x, p). Then £„(x), 
in Eq. ([/]), is the fluctuation of the current density 



Ut,x) = e/-^^£(A<-An«). 



3 The distribution function of anti-particles deviates from equilibrium by the amount — Sn(t, x, p). 
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One property of follows from charge conservation, d v £" = 0. To find the correlation 
of £ v (x), we first observe that the correlation of the thermal fluctuations An±(i, x, p) at 
coinciding times is given by 

(A<(t,x,p)A<'(t,x',p')> = (An!(f,x,p)A/(t J x' I p')) = 

(2nf5(p - p')(5(x - x')^'n(|p|)(l - n(|p|)) (9) 

(the correlator (An+An_) vanishes). The spatial delta function <5(x — x') reflects the 
short- distance nature of the fluctuations included in An±(t, x, p). Since we are interested 
in the correlation of £ u (x) at different times, we also have to know the time evolution of 
An±(t, x, p). It is obtained by solving the following transport equationF] 

(! +v £) An(( - x - p)=o - < io) 

Its solution is An(t, x, p) = An(0, x — v£, p). From this we infer that the noise correlator 



is 



*)Ut', y)> = 4e 2 JJ^ v»vj{x - y - v(t - t'))n(\ V \){l - n(|p|)) . (11) 

The chain of arguments which lead to Eq. ( |XT|) is somewhat intuitive. We can, however, 
verify this result in various ways. One way is to compute the correlator directly in terms 
of the fundamental fermionic field operator, if), 

(U x )Uv)) = e 2 mx)j^(x) ■ i>(yh„i>(y)) (12) 

The RHS of Eq. (|l^) is easily evaluated in the momentum representation. It is, in the 
approximation of free fermions, 

(£>, qje.K q')> = 4e 2 / ^ v,v v 5{u - vq)n(|p|)(l - B(|p|)) 

x (2tt) 4 (5(cj - u')5{q - q') . (13) 
It is simple to show that Eq. flT3] ) is just the Fourier transform of Eq. ([II]). 

4 This equation is appropriate for describing the propagation of fluctuations on a time scale much 
smaller than the transport mean free time. 
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We can also verify Eqs. ( [TT|JTB| ) using the fluctuation-dissipation theorem. Schemat- 
ically, we can solve Eq. (|6|]7]) to relate the power spectrum of the transverse part of 
with that of £ M by the very simple formula, 

( a>, W ^)>- <^;g (14) 



We did check that Eq. ([14]) leads to the spectrum of ^-fluctuations satisfying the 
fluctuation-dissipation theorem, Eq. (ffl). 

Having obtained the correlator, we now proceed to developing an effective dynamics 
for the long range modes of hot QED. We can solve the Vlasov equation Eq. (^) for 
<5n(£,x, p), and use it to obtain an equation for the induced current responding to the 
external field A^. We obtain 

5ra(t,x,p) = -e TduvE(t-u,x-uv) dn ^ P ^ , (15) 

n "\ PI 



anc 



x) = 4e 1 7% v5n(t ' x ' p) = ~tt~ / ^r^l^ - y) • E (* - l x - yl. y)l • ( 16 ) 

J (27t) ,:1 127T j |x — y| 4 

Until now, our only assumption has been that to and k are much smaller than T. In 
the regime to ~ g 4 T, k ~ g 2 T, additional simplification applies. In Eq. (|16|) we can replace 
t— |x — y| with t, since the typical time scale is (g 4 T) _1 and x — y ~ (g 2 T) _1 <C ((yf 4 T) _1 . 
Substituting (|IED into Eq. (|7|) and recalling that E, = A, is negligible in comparison to 
V x B, where B = V x A is the magnetic field, we obtain an equation which evolves long 
range fields in hot QED, 



V x [V x A(t, x)] = ^ / dy £-lL [(x - y) • A(f , y)] + £(t, x) . (17) 



Finally, we observe from Eq. (|Tl|) that the noise correlator is only non-vanishing on the 
light cone, t = ±|x|. Since x is typically much smaller(by a factor of g 2 ) than the time 
scale we are interested in, we can approximate the time dependence of the noise correlator 
with a delta function. To find the spatial dependence of £„, we integrate the RHS of Eq. 
(|TT| ) with respect to t and find, 

<6(f,x)&(f,y)> = ^ {x ~ yUx ~ y)j 6(t-t'). (18) 
07T |x — y| 4 



'Since the important modes are magnetic ones, we can limit ourselves to the three-vector j. 
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Eqs. ( p] , lI%D are the final equations describing the evolution of the soft modes in hot QED. 

Before going further, we need to justify the neglect of the collision term in the Vlasov 
equation. The typical time scale for scattering to significantly affect particle distribution 
functions, is determined by the transport cross section and is of order (g^T)^ 1 . Therefore 
one would think that a scattering term should be included in the Vlasov equation (Eq. (j^)) 
since this time scale is of the same order as the time scale under investigation. However, 
this is not so as we now explain. According to Eq. (|T6|) , the current density depends 
on the background field A^{y) at time moments yo so that xo — yo ~ (g 2 T)~ 1 . This 
means that, on the time scale of (g 4 T) _1 , the plasma responds almost instantaneously to 
a change in the soft background, and the slow relaxation due to particle scattering has 
only a negligible effect. We can verify this statement by including a simple collision term 
to the Vlasov equation, 

at + v 9xj Mt ' x ' p) + ^ Ev ^f = ~r 5n{t ' x ' P) ' 

where r is the transport mean free time ~ (g 4 T) _1 , and observe that its sole effect is 
the smearing of the delta-function over u — vq ~ g 4 T in Eq. (|I^). This has a negligible 
effect on the final equation (|17|) and the noise correlator (|18|) (the effect becomes large 
only when considering much smaller spatial momentum scale, g 4 T). As a result, we can 
neglect the effect of scattering on the behavior of the plasma since the scale of spatial 
momentum responsible for non-perturbative physics in hot QCD has been shown to be 

irT. 

4. By now, we have accumulated all the tools to finally address the dynamics of 
the soft modes in a non-Abelian plasma in the non-perturbative regime of interest. We 
will consider an SU(iV) pure gauge theory. The distribution function of gluons in a 
plasma is a matrix N ab (t, x, p) in color spacef]. The mean value of N ab is the equilibrium 
distribution ra(|p|)<5 af> , and, as before, fluctuations of N ab can be divided into soft and 
hard components. The soft component can be written in the form —if abc Sn c where Sn is 
the parameter entering the Vlasov equations (|J ||) . The hard component AN ab gives rise 
to the stochastic source in the RHS of Eq. (|5|), 

= ^/^^r bc E^ c (t,x, P ) (19) 



a 



denotes gluon polarization). The one-time correlator of AN ab is given by a trivial 



D We assume that the distribution is trivial with respect to gluon polarization, see, for example, 
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generalization of Eq. (^) , 

(AN?(t, x, p)A</(t, x', p')> = (2n) 3 5(p - p')5(x - x')rt 6b V Q 'n(|p|)(l + n(|p|)) . 

(20) 

Furthermore, AN ab satisfies the transport equation, 

v^V^AN = (21) 

where = + [A^, . . .], .4jf = gf abc A c ^. The substitution of the covariant derivative 
to the space derivative d x in Eq. (0) reflects the fact that hard gluon precesses in color 
space when propagating on the soft background. Remembering that the field strength A 
is typically of order gT, we observe that d x is of the same order as gA ~ g 2 T, in which 
case, the gA in the covariant derivative term cannot be dropped.^] Combining Eqs. (19). 
and (pip, we find for the noise correlator, 



(^,x)^(t',y)) = 2/iV^(t,x;t',y) |^^(x-y- v(f-f ))n(|p|) (1 + n(|p|)) . 

(22) 

U (x, y) is the Wilson line connecting x and y in the adjoint representation 

U(x,y) = Texpl- fdz^An ( z ) 
V y 

The presence of the Wilson line renders Eq. (0) explicitly gauge covariant. 

The generalization of Eqs. (0,0) is straightforward. First, we solve the Vlasov equa- 
tion Eq. (fD with respect to 5n a (t, x, p), 

5n a (t, x, p) = -ff ^pp Jdu ^ a6 (^, x - w) v • E b (x - w) . (23) 
IPI 

We then insert the expression obtained in the expression for the induced current which 
appears on the RHS of the field equation Eq. (H), 

#(f , x) = 2gN J j^- 3 vM a (t, x, p) . (24) 



7 In more physical terms, the color orientation of a hard particle changes essentially when moving over 
distances of order (<7 2 T) _1 . In contrast, the deviation of the particle trajectory from the straight line 
on the scale of (g 2 ! 1 )" 1 is negligible. In fact, the magnetic field B ~ dA ~ g 3 T 2 corresponds to the 
curvature radius of the trajectory of hard particles of order T/gB ~ (g 4 T) _1 which is much larger than 
the scale we arc interested in. 
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The equation obtained is non-local in space and time. However, we have already argued 
that, in the regime ui ~ g A T, k ~ g 2 T, the non-locality in time is confined to a time 
interval of order ~ (g 2 T)~ 1 , much smaller than the typical time scale of variation of 
the field ~ (g^T)^ 1 . Hence, in the limit 5 < 1, we can ignore the non- locality in time. 
Inserting Eq. fl24j) in the field equation (Eq. (|)), and using the fact that E can be neglected 
in comparison to D x A, we find, after little algebra, 



a 2 T 2 N r x — v 

[D x [D x A(f , x)]f = J dy —-j^ U ab (t, x, y) A% y) ■ (x - y) + £ a (t, x) (25) 



with the noise correlator ( p2|) taking the form 



m,x)£(t'iy)) = 9 -^U*%x,y) (x y)i{x u y)j 6(t-t'). (26) 



Gtt " |x — y 



It can be checked that Eq. fl2"6|) is consistent with the transversality condition = 0.0 
Equations (^) and (|26"D , constitute the highlight of our paper. 

5. In summary, we have derived a Langevin-type equation, Eq. (|25|), which describes 
the non-perturbative behavior of very soft modes in hot gauge theories. They are non- 
local in space, but local in time. Equation ( p5|) is a first-order differential equation with 
respect to t, and A appears in the equation in a linear fashion (there is no terms like A 2 ). 
The noise correlator is local (white) in time, but has a non-trivial spatial dependence. Eq. 



(25) can be used for simulating the dynamics of soft modes in a variety of contexts. To 
do so, one generates, at any given time, a noise £, transverse and distributed according 
to Eq. (^), then runs the system to the next time-slice using Eq. fl25|) ; one repeats these 
steps until the system reaches a stationary state. The benefit of our formalism is that 
one can concentrate on relevant modes of spatial momenta of order g 2 T only, since all 
higher modes have already been integrated out in deriving Eqs. (25], 26). Therefore, one 
can use a lattice with the lattice cutoff a -1 of order g 2 T, but not T or gT as in most of 
the methods suggested so far (the hard modes have been taken into account analytically), 
which leads to the possibility of doing calculations with a much larger lattice spacing. 
However, to actually use Eq. (^) for simulating the dynamics of soft modes, one has 
to invert the linear operator acting on A on the RHS of Eq. (P5[). This seems to be a 
potentially non-trivial problem for numerical calculations, and we hope that this obstacle 
can be overcome. 



This condition is required for Eq. (E5|) to have solutions with respect to A. 
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Application I: Baryon number violation at high T 



1. Our main motivation for developing an effective theory for soft modes is their involve- 
ment in numerous important phenomena taking place in a high temperature environment, 
one of which being the mechanism of baryon number violation at high T. Existing numer- 
ical simulations, and in particular the most recent ones [l], |), suggest that, T w , the rate 
of baryon number violation per unit volume, is ko*T with k ~ 1. This corroborates the 
most naive dimensional estimates which suggest that T scales as ~ £ 4 where £, ~ (g 2 T)~ 1 , 
is a natural length scale of the problem. Two reasons suggest caution: 

- 1 - First, the above naive dimensional analysis does not recognize that the relevant time 

scale is (g 4 T) -1 rather than (g^T) -1 . 

- 2 - Second, existing numerical simulations incorporate high momenta modes classically. 

This method is potentially sensitive to an ultraviolet cutoff and needs to be imple- 



mented with great care in order not to induce lattice artifacts) 10 . 



The effective theory developed in the previous section, has the virtue of incorporating 
the physics of the hard modes into the physics of the soft modes, hence, it alleviates the 
sensitivity on very short distance fluctuations. The outcome is a Langevin-type equation 
(one derivative in time, two derivatives in space) which suggests a cross-barrier transition 
of a diffusive type with natural length scale and time scale of order (g 2 T)~ 1 and (g 4 ! 1 ) -1 , 
respectively. 

We are confident that our results will become a tool to further study the rate of baryon 
number violation at high T. The qualitative and quantitative understanding of the latter 
has important consequences, some of which are illustrated below. 

2. Let us consider the outset of electroweak baryogenesis. The most recent compu- 
tation of the baryon asymmetry produced at the electroweak era, in the framework of 



some super symmetric theories 0], concluded that in a region of parameter space where 
all squarks are heavy and nearly degenerate (> 100 GeV), the baryon asymmetry scales 
as sin 8 cp T w /T s . Here, T s is the rate per unit volume of violation of the axial charge 
due to strong anomalous QCD processes and 6 cp is a CP violating phase. In contrast, 
in a regime of non- degeneracy, with, say, the super-partners of the third generation of 
quarks lighter than T and lighter than the other squarks, the baryon asymmetry scales 
as sm9 cp T w . 

In the speculative situation where T WjS scales as s rather than as a 4 s , the result 
of these computations is offset in respect to those quoted in Ref. WM, by a factor of 



ctw/cts ~ 1/3 in the degenerate case and a factor of a w ~ 1/30, in the non-degenerate 
case. In the latter case, these corrections are very significant as they correspondingly shift 
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the quoted lower bounds on the CP violating phase 9 cp from a range 10 -4 — 1CT 2 to a 
range ~ 30 larger, bringing these bounds to a "dangerous" overlap with the quoted lower 
bounds obtained from electric dipole moment experiments ||15|| . 

3. Finally, let us briefly speculate on another possible interesting consequence for 
electroweak baryogenesis. Another source of significant constraint on electroweak baryo- 
genesis is the preservation of the baryon asymmetry at times subsequent to its production. 
This requires that baryon number violating processes be suppressed in the broken phase 
after the transition. Currently, this requirement yields very strong constraints on the pa- 
rameters of the Higgs sector. It remains to be explored whether the above considerations 
will affect the determination of these constraints. 



Application II: Parameters of the electroweak phase 
transition 

Let us now turn to a second important application of our formulation of long range 
physics in a non-Abelian plasma: the determination of the parameters of the electroweak 
phase transition. We refer, more specifically, to the order of the transition, and, if first 
order, to the nucleation rate, and the dynamical properties of the phase interface such 
as final Higgs expectation value, shape, profile and velocity. The knowledge of these 
quantities is required for a quantitative study of electroweak baryogenesis. One of the 
computational tools often used in their calculation is the equilibrium free energy, V(4>, T), 
of the system. The perturbative determination of the latter is known to be plagued with 
infrared divergences. Those arising from the physics at scale gT can be accounted for 
perturbatively by resummation of hard thermal loops 1 16, 17]. Others, arising at scale g 2 T 



in the magnetic components of the gauge sector, remain unaccounted for perturbatively; 
their presence is an obstacle to a complete understanding of the outcome of the electroweak 
phase transition. 

The stochastic dynamics of soft modes presented above sheds new light on these is- 
sues and provides a new tool for resolving them. To illustrate how, we begin by briefly 
presenting a real-time calculation of the free energy V((f),T). 

1. We consider a gas of gauge bosons with mass m = g<p/2, in the presence of a space- 
dependent Higgs background (f>(x) which interpolates smoothly between phases of broken 
and unbroken electroweak symmetry: 0(— oo) = 0, <j)(+oo) = <po and <9 x </>(±oo) = 0. 
Focusing on the gauge components of the plasma and their interactions with the kink, we 
can simultaneously describe the dynamics of the kink and of the particles distributions 
n a (x, p) with the following equations |l8f 
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d a , \ 1 dm 2 (x) d n . . . 



,2/ 



n , x dV(d>) dm 2 (x) r do 1 „. . . 

^ x) = --M + -±l J^L-^j,), (28) 

where e p is the space-dependent energy of a particle with momentum p at point x. 
The force term in the first Vlasov equation above is derived from energy conservation, 
which, in the kink frame, implies p 2 + m 2 (x) = constant. We assume the existence of a 
stationary situation. This is guaranteed dynamically by the presence of the collision term 
C(t). Without the latter, the kink would continuously accelerate leading to an unpleasant 



situation 18, 17 



The above equations supplemented with appropriate boundary conditions for the kink 
and with e p (— oo) = p 2 , not only provide a unique solution to the propagation of a 
phase interface (in a regime where the curvature can be ignored) such as its velocity 
relative to the plasma^ and its profile, but also provide us with a powerful scheme for 
computing V(<p,T). A useful quantity to look at, is the space variation of T tlv (x), the 
stress energy tensor of the system. The leading order term of its expansion in Vk, the 
relative kink-plasma velocity, yields the free energy of the system T xx (x) — T xx (—oo) = 
V(0,T)-V((f)(x),T)+v k L((j)(x),T) + O(v 2 k ). Because the "drag force" L{<f>{x), T) vanishes 
as the collision term in the first Vlasov equation vanishes, the effective potential is the 
solution of Eq. (p8|) in the collisionless limit. The solution is £p(x) = p 2 +m 2 (x) and, 
V(<f>,T) given by the so-called "one-loop" effective potential 0]. Its minimum, O , scales 
as 0o ~ gT m^/m^ which, in turn, fixes the size of the kink to be of the order of 
Lk ~ / dtp/ JV{4>, T) ~ 1/(? 2 T {m H /m w ). This approximation of the free energy has 
been known for a long time to be inadequate as it does not incorporate any of the long 
range plasma physics. In particular, in the limit (f) — > 0, there is an uncontrollable number 
of soft quanta ~ T/g<p whose interactions (with strength g 2 ) among themselves and with 
the hard modes yield uncontrollably large corrections to the free energy. The remaining 
of this section addresses that question. 

2. Inclusion of the gT -scale physics. Plasma physics at the scale ~ gT ("plasmon 
physics") becomes relevant when g<p ~ gT, that is for <ft < T. From our estimate of </>o 
above, we see that this physics affects the determination of the free energy for all values of 
< 00; i-e., across the whole kink. The resolution of these difficulties is know|16, 17]. One 



has to incorporate hard thermal loops in the calculation of V(0, T). In the above real time 
framework, this operation can be performed very simply by noting that plasmon physics is 



3 Eqs. (g8|) were used in Ref. to compute the relative motion of the kink in the plasma. 
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contained within a distance ~ 1/gT < ~ l/g 2 T: the plasmon physics is "local" inside 
the kink and is in a large extent independent on the properties of the latter. One can 
then describe the kink as coupled to a gas of free bosons with longitudinal modes replaced 
with plasmons. The dispersion relation of a plasmon is, to a sufficient approximation: 
£p = p 2 + Il(p = 0), where yTl(p = 0) is the plasmon frequency, ~ gT. Incorporating this 
modified dispersion relation in the boundary conditions at x — > — oo, yields an expression 
for the source s(x), which, after simple algebra, reproduces a form of the free energy known 
in the literature as the "one loop improved" effective potential |T7[. This is a significant 
improvement because gauge interactions among plasmons are computable perturbatively, 
their effective strength scales according to g 2 T / \Jll(p = 0) ~ g < 1. There remains an 
infrared behavior in the transverse components of the gauge sector. 

3. Inclusion of the g 2 T -scale physics. It is transparent from the above analysis that 
the calculation of V{4>,T) involves the natural scaled ~ l/g 2 T which is of the order 
of the one characteristic of the non-linearity in a hot gauge plasma. Any determination 
of V(4>,T) which claims to incorporate the latter physics, is doomed to be intertwined 
with the physics of the g 2 T modes, in contrast with the situation we encountered with 
the gT modes. A path for a resolution is provided by evolving simultaneously the plasma 
in a Higgs background and in a soft gauge background. The equations consist in a 
generalization of Eqs. (p5|) and (|26|) to the presence of a space-dependent background 
Higgs field, that is, a Vlasov equation supplemented with the gauge field equations and 
the Higgs equation of motion. In contrast with the situation with the gT-modes, these 
last two equations do not decouple from each other because the typical lengths of variation 
of both backgrounds are comparable, ~ (g 2 T) _1 . In particular, the space derivative on 
the LHS of Eq. (j28D is to be substituted with a covariant derivative. These equations, 
whose study is beyond the scope of this letter, constitute a natural real-time framework 
to incorporate the soft gauge physics in the calculation of the thermodynamics of the 
electroweak phase transition. 
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